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Abstract. 

It is well known that the Solar System is presently moving through a partially ionized local interstellar medium. 
This gives rise to a counter-flow situation requiring a consistent description of behaviour of the two fluids - ions and 
neutral atoms - which are dynamically coupled by mutual charge exchange processes. Solutions to this problem 
have been offered in the literature, all relying on the assumption that the proton fluid, even under evidently 
nonequilibrium conditions, can be expected to stay in a highly-relaxated distribution function given by mono- 
Maxwellians shifted by the local proton bulk velocity. Here we check the validity of this assumption, calculating 
on the basis of a Boltzmann-kinetic approach the actually occurring deviations. As we show, especially for low 
degrees of ionization, £ < 0.3, both the H-atoms and protons involved do generate in the heliospheric interface 
clearly pronounced deviations from shifted Maxwellians with asymmetrically shaped distribution functions giving 
rise to non-convective transport processes and heat conduction flows. Also in the inner heliosheath region and in 
the heliotail deviations of the proton distribution from the hydrodynamic one must be expected. This sheds new 
light on the correctness of current calculations of H-atom distribution functions prevailing in the inner heliosphere 
and also of the Lyman-a absorption features in stellar spectra due to the presence of the hydrogen wall atoms. 
Deviations from LTE-functions would be even more pronounced in magnetic interfaces, which via CGL-effects 
cause temperature anisotropies to arise. 



1. Introduction 

The problem of the heliospheric interface, where H-atoms 
and protons are effectively coupled by charge exchange in- 
teractions, has often been faced in the literature. In gen- 
eral, it was recognized very early that the passage of neu- 
tral interstellar atoms (O, H) through the plasma inter- 
face ahead of the solar wind termination shock needs a ki- 
netic treatment, since the relevant charge exchange mean 
free paths between H-atoms and protons are comparable 
to or even larger than the typical structure scales of the 
interface plasma flow (i.e., Knudsen numbers are equal 
to or smaller than 1; see Ripken & Fahr 1983; Fahr & 
Ripken 1984; Fahr 1991; Osterbart & Fahr 1992; Baranov 
& Malama 1993 ; McNutt et al., 1998, Bzowski et al. 2000; 
Izmodenov et al. 2001) . Nevertheless, due to mathemati- 
cal complications associated with such kinetic treatments 
of the problem, many heliospheric models have appeared 
in the literature which use hydrodynamic treatments of 
the two fluids, H-atoms and protons, coupled by charge ex- 
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change reactions (for recent reviews, see Zank 1999; Fahr 
2003b). 

The hydrodynamics applied in all these approaches is 
restricted to the description of the space-time behavior of 
the lowest hydrodynamic moments like density, bulk ve- 
locity and pressure and the local distribution function is 
taken to be solely a function of these moments in the form 
of shifted Maxwellians. As suggested in Fig^ the perma- 
nent supply of newly charge- exchanged particles into the 
local distribution functions will maintain the resulting dis- 
tribution (shown by dashed lines) far from a three-moment 
HD distribution. 

This has been clearly recognized by Baranov and 
Malama (1993, 1995), who for this reasons perfected a 
kinetic treatment for the neutral H-atoms, even in view 
of the mathematical complications. The semi-kinetic ap- 
proach offered by them, though treating the H-atom ki- 
netically, is still based on the assumption that the protons 
can be described as a hydrodynamic fluid. But this is not 
true under special conditions, as we will demonstrate here. 
The protons generate deviations from a 3-moment HD- 
distribution. Therefore to represent the two fluids interact- 
ing by charge exchange reactions one would need a kinetic 
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Fig. 1. Schematic view of the "old" proton and H-atom 
distribution functions (fo lP and /o,ff, light gray and dark 
gray, respectively), with different bulk velocities u p , uh 
and temperatures, interacting with each other by charge 
exchange. The charge exchange gives rise to secondary 
populations of protons and H-atoms that inherit the kine- 
matic properties of their respective source populations 
of unperturbed H-atoms and protons. As a result, non- 
Maxwellian distribution functions f p , fn of both protons 
and H-atoms appear, shown in black lines. 

description both for the H-atoms and the protons. This 
highly complicated model will not be offered in this paper, 
but we present calculations which clearly make visible the 
resulting deviations of H-atom and proton distributions 
from 3-moment HD distributions. One other point, which 
was already emphasized by Fahr (2003a, b), concerns the 
fact that hydrodynamic two-fluid descriptions of the in- 
terface flows use charge-exchange coupling terms which 
are only justified in cases of supersonic bulk velocity dif- 
ferences. Under realistic interface conditions the resulting 
bulk velocity differences, however, appear to have subsonic 
magnitudes and thus require revised formulations of the 
charge exchange coupling terms. In the following we thus 
subject hydrodynamic approaches to a kinetic control. 

2. Theoretical approach 

We assume that protons and hydrogen atoms in the in- 
terface are coupled solely by charge exchange processes. 



Neutral density 



0.12 




0.1 I 1 1 1 1 1 J 

120 140 160 ISO 200 220 240 
r[AU] 
Neutral bulk velocity 





-10 




-12.5 




-15 


1 


-17.5 








-20 




-22.5 




-25 



120 140 160 180 200 220 240 
r [AU] 

Neutral thermal velocity 




11.5 



120 140 160 180 200 220 240 
r[AU] 

Collisional relaxation time 



7xl0 9 




120 140 160 180 200 220 240 



r[AU] 



Fig. 2. Macroscopic parameters of H-atoms between the 
bow shock and the heliopause, calculated with the use of 
the multi-fluid hydrodynamical treatment developed by 
Fahr et al. (2000). The bow shock is at - 250 AU, the 
heliopause at ~ 125 AU. 
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Fig. 3. Macroscopic parameters of protons between the 
bow shock and the heliopause, calculated with the use of 
the multi-fluid hydrodynamical treatment developed by 
Fahr et al. (2000). 



For simplicity, we limit ourselves to calculations on the 
stagnation line only. Denoting the hydrogen distribution 
function by fjj and the proton distribution function by 
f p we furthermore assume that the two distribution func- 
tions can be represented by Maxwellian cores fo.n, fo,p 
from hydrodynamic calculations and by non-relaxed cor- 
rections fi t H, fi,p due to mutual implantations of new 
charge exchanged particles, as shown in the following for- 
mulae and qualitatively in Fig^ 



f H (v, s) = fo tH (v, s) + f x ,H (v, s) 
f P (v, s) = f , p (v, s) + h.p (v, s) 



(1) 
(2) 



where s measures the spatial position on the symmetry 
axis. At this symmetry axis, and for stationary conditions, 
we have the following system of Boltzmann-kinetic equa- 
tions in phase space, which couple the protons and hy- 
drogen atoms via mutual losses and gains (see, e.g., Fahr 
1996): 



v cost 



ds 



S + fH 5 f H foM - Ih 



St 



5t 



TH 



dfp _ S + f p 8 f p f ,p - f P 



ds 



St 



St 



(3) 



(4) 



In these equations s is a line element along the axis in the 
interface and tjj, t p are relaxation times, separate for pro- 
tons and hydrogen atoms and dependent on s and the par- 
ticle velocities. The terms with 8 + and S~ are, correspond- 
ingly, the production and loss terms, as shown in formulae 
13 and |SJ and the remaining terms in Eqs |3] and 0] are the 
relaxation terms operating to restitute the hydrodynamic 
solutions for appropriate distribution functions. The actu- 
ally resulting distribution function / which, had the local 
charge exchange influences been stopped, would relax to- 
wards the hydrodynamic core distribution /o within one 
relaxation period. But as long as charge exchange is oper- 
ating, the resulting distribution is permanently kept off 
this function Jq. In other words, we assume that two- 
fluid hydrodynamics would deliver correct results in the 
interface had the H-atom and proton relaxation processes 
operated fast enough, i.e., if th , t p << r ex , r ex being the 
charge exchange period between the H-atoms and protons. 
Only because th, t p ~ r ex deviations from the hydrody- 
namic core functions are pronounced. The above system 
of differential equations is valid for an arbitrary velocity 
vector v (v, cos 9) but for a proton velocity vector identical 
to the H atom velocity vector. 9 is the angle between the 
direction of atom's motion and the inflow axis z (and con- 
sequently v cos 9 = v z ). No matter what the forms of the 
distribution functions are, a phase-space hydrogen gain 
is then a phase-space proton loss and vice versa. Hence 
equations |3] and 0] are modified to 

dfn _ S~f p S~f H fx,H 



v cos 9 



ds 



St 



St 



a d fp 
vcos9 —f- = 

ds 



S fn _ S fp 
St St 



TH 

fl,p 



(5) 



(G) 
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Specifically, the two remaining coupling terms are of 
the following forms: 



= f H (v,s,e H )x (7) 
x I v re i (v,v p ,9 p ) f p (s,v P7 9p) a (v re i) d 3 v p 



= fp(v,s,6 p ) x (8) 
x v re i(v,v H ,9 H ) f H (s,v H ,9 H ) a(v rel ) d 3 v H 



S-f H 
St 



s-fp 

St 



Here v re i (v,v p ,9 p ) denotes the relative velocity between 
an H-atom of velocity v and a proton of velocity v p inclined 
at an angle 9 P to the proton bulk speed u p . a (v re i) is 
the H-atom - proton charge exchange cross section for 
particles with a relative velocity v re i . 

The core functions fo,H/p m the integrands will be 
taken from a hydrodynamic multifluid simulation code 
(Fahr et al. 2000) describing the interface system. They 
are functions of s and defined as follows 



fo,H {v,s,9 H ) = 



1 



7T C H 



x exp 



n H x 



H 



(9) 



2uh v cos 9h 



fo, P (v, s, Op) 



7T Cr 



x exp 



H 



u p — 2u p v cos ( 



(10) 



where uh, u p are, respectively, the local bulk velocities of 
hydrogen atoms and protons, which depend on the posi- 
tion and thus here are given as functions of the path length 
s measured between the outer shock and the heliopause; 
nu, n p are number densities, also dependent on s; 9h, 
9 P are angles between the local bulk velocities, which of 
course on the axis must be parallel to the inflow axis, and 
the individual velocity V] and ch , c p are thermal velocities 
of the core hydrogen atoms and protons, also dependent 
on s via temperatures Th , T p : 



a _ 2K T H (s) 

C H — 



rn 



2 Ph (f) 
m uh (s) 



2K T p (a) = 2 P p (s) 
m m n p (s) 



(11) 



(12) 



We assume that the proton mass m is equal to the hy- 
drogen mass. The temperatures are also functions of s. In 
Figures |21 and 13 we present the core hydrodynamic quanti- 
ties (density, bulk velocity, thermal velocity, and tempera- 
ture) as functions of s, as they result from the multi- fluid 
interaction code developed by Fahr et al. (2000), which 
here was run for the following input parameters: solar wind 
speed: 400 km/s, solar wind density at 1 AU: 5 cm -3 , solar 



wind temperature at 1 AU: 10 5 K, interstellar gas speed: 
26 km/s, interstellar proton density: 0.1 cm~ 3 , interstel- 
lar hydrogen atom density: 0.1 cm~ 3 , interstellar gas tem- 
perature: 8000 K. The influences of pick-up ions and of 
anomalous and regular components of cosmic rays were 
switched off in this modelling for consistency with our ki- 
netic approach presented in this paper, where couplings 
to these particle populations are not taken into account. 

We thus start from a two- fluid hydrodynamic interface 
model in which H-atoms are represented by one single fluid 
only. This is different from modellings like that presented 
by Zank et al. (1996), where neutral hydrogen is decribed 
by three different fluids representing H-atoms a) coming 
directly from the interstellar medium, b) originating in the 
hclioshcath, and c) originating in the inner heliosphcrc. 
For our approach here, where we want to check on kinetic 
deviations from hydrodynamic H-atom fluid approxima- 
tions, we consider it to be more practical to start from a 
simple mono-Maxwellian representation of the H-atoms, 
thereby making our concept logically more convincing. 

We rewrite the coupling terms defined in equations 
and[S]so that the core + correction form of the distribution 
function (see Eq^and^J is taken into account: 

" Jn f H (v,s,e H )x (13) 



St 



+ 



V re l (V, V p ) / ,p (s, V p ) G (v rei ) d 3 Vp + 



Vrei {v, V p ) fi >p (s, v p ) a (v re i) d?v p 



s-f P 

St 



fp(v,s,9 p ) x 



(14) 



X / Vrel (V,V H ) fo,H {s ,V H ) <J {v r el) <P V H + 



el (V, V H ) fi >H (S, V H ) a (Vrel) d 3 V H 



and then we neglect the correction terms in the integrals, 
since they only contribute according to their partial den- 
sities with ni t H -C no t H, n\ tP <C no tP . We assume that 
the total density ri\ given by /i is negligible with respect 
to the total density no given by fg. This does not mean 
that fx itself is everywhere small with respect to /o, as 
shown in Fig^ where velocity regions where f\ > fo can 
be seen. We thus assume that in Eas.HSIandlTBIonlv terms 
(/o,i + fi,i) J fo.j--- must be taken into account. Hence 
only the shifted Maxwellians stay in the integrands of 
equations 1131 and 1141 yielding: 



St 



s-f P 

St 



= f„(v,s,6 H )x (15) 

X V re l(v,Vp) fo,p(s,V p ) O (v re l) d 3 V p 



= f p (v,s,9p) x (16) 

X V re l(v,V H ) f()M(s,V H ) <T (v re l) d 3 V H 
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The validity of this assumption can be checked apos- 
teriori in the results presented in Figures through [SJ It 
can be seen that close to the heliopause this assumption 
is violated. 

The dependence of the charge exchange cross section 
on the relative velocity of colliding particles can be de- 
scribed by the well-known formula given by Maher & 
Tinsley (1977): 



v{Vrel) = (A + B \og(v rel )) 2 



(17) 



where A = 1.6 10~ 7 and B = —6.8 10~ 9 . Now we further 
simplify, by assuming that the charge exchange cross sec- 
tion can be expanded into a constant + linear expansion 
term, i.e., expanding around the bulk velocity u p of the 
protons and bulk velocity Uh of the H-atoms, respectively: 



cr p (v re i yP ) = (To, p + <Ji,p {v re i lV — Vq) 



&H (v re l,H) = &0M + &1,H (v re l,H ~ Vq) 



(18) 
(19) 



We take as v appropriate bulk velocities u p , uh and we 
change variables so that now we will be in a reference 
system co-moving with the protons or vice versa with the 
hydrogen atoms. The coefficients Co and <j\ in equations 
1181 and [Tulare calculated from the following formulae: 





= (A + Blogu H ) 2 




(20) 


&1,H 


do . . „ , , 
= —(u H ) = 2(A 
av 


-\-Blogu H ) 

UH 


(21) 




= (A + B\ogu p ) 2 




(22) 


01,p 


dcr . . „ , . 
= ^ (Up) = 2(A + 


B log Up) — 

Up 


(23) 



where u p , ur are velocities normalized by the velocity 
of 1 cm/s. The relative velocity v re i in equations 1181 and 
1191 is expressed as v re i tP — \u p — v\ for protons and as 
v T el,H = \uh — v\ for hydrogen atoms. 

In the co-moving reference systems the Maxwellian 
core distribution functions are the following: 



f°* ( v ' s ) = ~, 7TT372 ex P 

(7T Cp (s)) 



c p (sY 



fo,H(v',s) 



n H (s) 



(7r c H (s)) 3 ^ 2 



exp 



ch (sY 



(24) 



(25) 



where v' are specific velocities in the appropriate co- 
moving frames. Evaluating equations 1151 and 1161 with the 
use of El El brings the following expressions for the cou- 
pling terms: 



St 



= fH 



0Q,p / VreLp fo,p d V p 



(26) 



+ 01,p / "rei.p U,p d 3 V p 



s-f P 

St 



fp <T0,H J V re l t H fo.H d 3 v' H - 

a\,H J v 2 relH f ,H d 3 v' H 



(27) 



The relative velocities under the integrands are defined as 
follows: 



^rel,p 



V 



rel,H 



— v + u H — 2v uh cos ( 



(28) 
(29) 



From Fahr & Mueller (1967) and Ripken & Fahr 
(1983), the mean relative velocities coming out from the 
first integral in equations l26l and l27l are given by formulae: 



{VreL; 



■. exp 



+ u v 1 + 



2u r 



erf 



(30) 



{Vrel,H) = —f= exp 
\/7T 



4 1 



H 



-UH 1 



For the second integral in formulae 1261 and 1271 we notice 
that the integration over 9 will bring to the terms with 
cos 9 and that - since u p , uh are local quantities which 
do not depend onv- the second term from the relative 
velocity expansion ( equations 1281 and 1291) will simply yield 
u 2 times the result of integration of the core distribution 
function over the velocity space. We therefore can denote 
this result simply by no h and n p , respectively. The re- 
maining integrations of the mean square velocity with re- 
spect to the core distribution functions yield simply the 
local thermal velocities and thus finally the integrations 
defined in equations [SB] and |23 lead to the following re- 
sults: 



St 



fn [(To 



01, 



!)] 



(32) 



S~ f 

— TT^ = U H f p [<J 0t H (VreLH) + 01, if ("t,H + U h)] ( 33 ) 

with the mean relative velocities defined in Equations 1301 
and|^ The set of differential equations, originally denned 
by formulac|2land21 is now obtained in the following form: 



v cos 9 



dfiM 
ds 



Q dfa,H fi.H ts . , s 

= -UCOS0— : n p (f ,H + Ji,h) x 

as th 

X [oo,p (Vrel.p) + 0"l,p (c^ + Mp)] + (34) 

+ n H (f . p + fx, p ) x 

X [a-Q.H ( V reLH) + Cri iH (c 2 H + U 2 H )] 



v cos 9 



df i, f 



= — v cos 9 



dfx,p fip / . » % 
- n H (fo, P + h, P ) x 



ds ds 

X [(Ta.H (Vrel,H) + a \,H ( c ff + u h)\ + ( 35 ) 

+ n p (f 0tH + f 1<H ) x 

x [(To. p (v re i iP ) + oi.p (cp + u 2 )] 
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The above equations can be rewritten so that the terms and analogously for the function /o, P - The terms of equa- 
multiplied by f\, p -H are collected together: tion0B|are evaluated as follows: 

dfo,H _ fo,H 
-i' cost) — (36) driH 



ds 



ds T H 

fi,Hn p (cr 0:P (v re i tP ) + ai tP (Cp + u 2 p )) + dfo.H 

fi, P n H [cr 0M (v TeljH ) + <j x ,h (c 2 H + u 2 H )] - du H 

fo,H n p [<T ,p (VreLp) + CTl.p (c 2 p + 1%)] + g j q h 

fo,p n H [°0M {Vrel,H) + °\,H (c H + u h)] 9c H ' C H 



fo.H 



u H 



H 



Ujj — 2v uh cos( 



H 



(47) 
(48) 

(49) 



v cos 9 



dfi.p 

ds 



= -vcosv — — — — (37) 

ds T p 

- fl, p riH (o"0,H (VreLH) + a l,H (<?H + u h)) + 
+ fl.H Tip [<T ,p (Vrel.p) + (T\,p (c p + U 2 )] - 

- fo.p n H [(TQ,H (VreLH) + C%ff {<?H + u h)] + 

+ fo,H n p [a , p (vrei.p) + <Ji,p (cp + Up)] 

These equations hence formally represent the following lin- 
ear system of differential equations: 

//,V // A (s)+B(s) h.„ + C(s) h. p (38) 



V cos 9- 



ds 



VCOS 0^ = D(s)+E(s) fi, H + F(s) f h 



ds 



(39) 

with the multiplicative factors A, B, C, D, E, F dehncd as 
follows: 

% - (40) 



A(s) 



-v cos i 



-fo,H n p [<7 0tP (v re i iP } + ai t p (cp + u 2 )] + 

+fo,p n H [<TQ,H (Vrel.fl) + &\,H {<?H + U h)] 



B(s) 



h Up {(XO p {V re l,p) 



ai,p {c 2 p + up)) 



C (s) = n H [a 0< H (Vrel,H) + &1,H (c 2 H + Uff) 



D(s) 



— V COS0 r^- + 

ds 

+fo,H rip [<T 0:P (VreLp) + <?l,p (c 2 + V, 2 )] - 

— fo,p riH (v re l,H) + P\,H (c H + U H )\ 



(41) 



(42) 



(43) 



E (s) = rip [<7 . p (v re i tP ) + <Ji.p (cp + u 2 )] 



(44) 



F{s) 



n H (CT 



H (VreLH) + <J\,H [<? H + U h)) ( 4 5) 



The derivatives of the unperturbed distribution func- 
tions given by Eauationl9landlTT)larc calculated as follows: 



dfo,H dfo,H dn H df 0H du H df Q ^ H dc H 



ds 



dn 



H 



ds 



du 



H 



ds 



dc 



H 



ds 



(46) 



Collecting them all together we obtain the following ex- 
pression for the derivatives of the unperturbed distribu- 
tion functions: 



dfo.H 

ds 



fo,H 



ds 



,2 



du 



H 



ds 



(50) 



u H — 2vuh cos 9 

CH 



dc H 
ds 



dfo, P 
ds 



fo.p 



ds 

2 



V COS 



v + Up — 2v Up cos t 




(51) 



Now we need the relaxation times, used in equa- 
tions [S] and For the relaxation of the neutral hydro- 
gen gas, the only mechanism is elastic H - H collisions. 
The corresponding relaxation time is numerically given by 
Brinkmann (1970) and Fahr (1996) in the following form: 



th 



6.62- 10 1U 
riH\fTn 



6.62-10 

riH c H 




;.5i • io 1 

n H Ch 



(52) 



for nn expressed in cm -3 and ch in cm/s. The period th 
turns out to be very long in comparison to charge exchange 
periods r ex , i.e., th/t cx ~ 10. The relaxation time for 
protons is more difficult to calculate. If, however, the sec- 
ondary protons have a bulk velocity small enough to not 
be shifted outside the broad thermal envelope of the "hy- 
drodynamic" proton population (in other words, the net 
distribution function does not feature two distinct peaks), 
then the wave-particle interaction by a two-stream insta- 
bility should not be significant. If, in addition, the pre- 
existing turbulence level in the unperturbed LIC is small 
at the distance scale comparable to and smaller than the 
size of the interface, then also non-linear interaction with 
MHD waves (mainly Alfven waves) may be negligible and 
the only relaxation mechanism will be Coulomb collisions 
between the old and newly-created protons. A more in- 
depth discussion will be presented in Sec. 4 at the end of 
this paper. 

The proton relaxation time is defined as the time be- 
tween Coulomb collisions of a proton traveling at v and 
protons belonging to the core population, described by the 
proton distribution function fo being a Maxwellian shifted 
by Up with thermal velocity c p , and is obtained from the 
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mean free path with respect to Coulomb collisions Ac in 
the form: 



Ar 



1 



\n p a c (v re l,p) 



(53) 



where | is the energy transfer efficiency term, 
equal to the reduced mass of the colliding particles: 
2 m p m p / (m p + m p ) 2 . 

The angle-averaged mean momentum-transfer 
Coulomb collision cross section ac in Eq. 1531 is given by 
the textbook formula (e.g., Oraevskij 1989): 



G c = °c,_lA 

where A is the Coulomb logarithm defined as 

A = ln^ 
P± 

(Tc.± is defined by the formula: 
ocx = n P±. 

and p± in Eg 1561 is defined by 

e 2 

— =KT p 
Px 

thus yielding 

e 2 



K T, 



v 



so consequently 



7re 



(kt p y 



(54) 



(55) 



(56) 



(57) 



(58) 



(59) 



In Easl54lthroughl59l K is Boltzmann constant; (v re i }P ) is 
the mean relative velocity of the test particle (the newly- 
created proton) with respect to the unperturbed proton 
population described by /o, P , defined in Eg. 1301 and e is 
the cgs elementary charge. Deep inside the interface, at 
200 AU from the Sun, the proton density n p = 0.18 cm" 3 , 
proton temperature T p = 2.3 10 4 K, and the relative speed 
(v r ei,p) = 9 10 5 cm/s. In EqE^lA/) is the Debye length, 
defined by the formula: 



A, 



KT n 



4-7T n p e 




(60) 



which for the above-mentioned conditions in the helio- 
spheric interface is equal to ~2460 cm, and the Coulomb 
logarithm A, defined in Eq.JSBl to —24.25. Thus the total 
Coulomb cross section ac, defined in Eo l54l is equal to 
~ 4 10 -13 cm 2 . Hence the proton relaxation time within 



the outer interface, 
order of a year: 

t v ~ 3.1 10 7 s 



at — 200 AU from the Sun, is of the 



(61) 



These times t p become even much larger in the inner inter- 
face and in the heliotail. These estimates may express the 



fact that relaxation processes are fairly unimportant here 
as well for the protons as especially for the H-atoms be- 
cause the relaxation periods are comparable to the mean 
passage times or structuring times 



uh. 



dn t 



10 s s 



(62) 



P ds 



In the final section of this paper we shall, however, come 
back to this point. The dependence of all the above ref- 
erenced quantities on the axial coordinates counted along 
the axis in the interface is shown in Figs 01 

3. Results and discussion 

We solved the system of equations 1381 1391 fsepararatelv for 
each v) using the hydrodynamic core distribution func- 
tions given in Easl9land llOI with the parameters obtained 
from hydrodynamic simulations by Fahr et al. (2000) and 
shown in Figures |21 and |3J In the approach presented here, 
where H-atoms are represented by a mono-Maxwellian 
plus kinetic deviations, we can only represent solutions 
of the distribution function for positive values of H-atom 
velocities. In order to present solutions for negative val- 
ues of H-atom velocities we would need to integrate Equs 
El and 01 in the opposite direction (negative increments of 
si) starting from an inner boundary value for the H-atom 
distribution function /o,jj (soj — ^h)- Since our hydrody- 
namic H-atom core functions vanish for negative velocity 
values, we cannot use any reasonable boundary condition 
to calculate the H-atom distribution functions for negative 
velocity values unless we admit additional H-atom fluids, 
as practised by Zank et al. (1996). The results for densities 
n Pt H, bulk velocities u Pi h, and temperatures T Pt jj have 
been calculated on the basis of the assumption that the 
underlying distribution functions f Pt H are three-moment 
functions, i.e., shifted Maxwellians. In contrast, in the ki- 
netic theory presented in this paper we have calculated 
those distribution functions actually occurring in the outer 
interface under the action of charge exchange processes 
between the locally colliding H-atoms and protons. These 
results we show in Figures|3through|Hl In these figures we 
have displayed at several interface positions on the sym- 
metry axis one-dimensional cuts through the kinetically 
resulting distribution functions f p ,H, showing both the un- 
perturbed hydrodynamical distribution functions /o(p,ff) 
and the kinetically perturbed, actually resulting distribu- 
tion functions / (PiH ) = fo( P ,H) + fi( P ,H)- 

First in Fig. we display w 2 -cuts through the distri- 
bution functions fu and f p resulting at various positions 
s = Si in the interface, when proton-proton relaxation pro- 
cesses are suppressed, but atom-atom relaxation processes 
are working. As one can clearly see, pronounced deviations 
of both resulting distribution functions f(p,H) from simple 
shifted Maxwellians, derived by hydrodynamical studies, 
can be recognized. These are even better identifiable in 
Fig- EI where we have displayed the ratios fi( P: H)/ fo(p,H)- 
Switching on the proton-proton Coulomb collision relax- 
ation, as developed in Eq. EH and as shown in Figs and 
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IH1 reduces these deviations in the case of proton distribu- 
tion functions to much milder degrees. Thus, under LISM 
conditions as used here in this paper, the hydrodynami- 
cally derived distribution functions for the interface pro- 
tons may be considered as good approximations, whereas 
the hydrodynamic approach to the H-atom distribution 
definitely is a poor approximation. Thus, due to rapid re- 
laxation processes under the conditions adopted in this 
paper deviations of protons from hydrodynamically de- 
rived Maxwellians can be considered as fairly small. 

In contrast, however, if interfaces for different LISM 
conditions need to be modelled, a kinetic revision of such 
proton distribution functions becomes necessary. For in- 
stance, this would be the case for lower LISM proton den- 
sities, because the proton-proton relaxation rate/period 
drastically goes down/up. For example, a LISM proton 
density of n PiOC — 0.01 cm~ 3 would increase the relax- 
ation rate by more than a factor of 10 with respect to 
the one used here, thus meaning in consequence that non- 
relaxated proton distribution functions must be expected. 
This would become important in the inner heliosheath and 
in the hcliotail. This all the more will be true if magnetic 
fields control the interface MHD dynamics and by CGL- 
effects (Chew et al. 1956) may lead to pronounced tem- 
perature anisotropics. 

As one can notice in Figs and |H1 the closer one ap- 
proaches the Sun starting from the outer bow shock, the 
more the resulting kinetic distribution function f( Pt H) de- 
viates from the unperturbed hydrodynamic distribution 
fo(p,H)y even reaching nonlinear amounts of such devia- 
tions. The specific form of these deviations thereby are 
not systematic with the coordinate s and could only per- 
haps be studied in more detail by the analysis of the 
higher moments of the distribution function /(p m, like 
the heat conduction flow or the temperature anisotropies. 
Qualitatively at least one can say - as it was already antic- 
ipated in the qualitative results sketched in Fig.0- that 
the solar/antisolar wings of the distributions f( p ,H) are 
enhanced with respect to the unperturbed ones fo( p .H) > 
wherever the local bulk velocities are larger than 

those U(j) t H) ■ Under these conditions heat conduction flows 
result from the distribution functions f( Pt H)y which for 
protons are directed in the direction of uh — u p and for 
H-atoms in the direction u p — uh- 

This fact has two consequences for the modelling of the 
H-atom presence in the inner and outer heliosphere. One 
is that the distribution function fn, which describes the 
injection of the LISM H-atoms into the inner heliosphere, 
is different from what was assumed up to now. This thus 
will influence the predictions made with respect to the 
H-atom phase-space distribution in the inner heliosphere 
and the related predictions on the resulting spectral and 
integrated Lyman-a resonance glow (see, e.g., Scherer et 
al. 1999). Another consequence may be that calculations 
of the absorption spectra caused by hydrogen atoms in 
the hydrogen wall ahead of the heliopause (see papers 
by Linsky & Wood 1996; Wood & Linsky 1997; Gayley 
et al. 1997; Izmodenov et al. 1999) may need revision. 



For instance, the asymmetric shapes of the resulting H- 
atom distribution function presented in this paper lead 
to asymmetric spectral absorption features at Lyman-a 
wavelengths which could remedy the problems in the cor- 
rect representation of the observed stellar spectra (see, 
e.g., Izmodenov et al. 1999). 

4. Outlook: Competing proton relaxation 
processes in heliospheric interface 

In the above calculations we have made use of the as- 
sumption that relaxation of the actually formed proton 
distribution function towards an associated equilibrium 
distribution / with S (/ /o ln/ d 3 v) /St — operates 
on time scales comparable to charge exchange time scales 
and the typical interface passage times. As proton-proton 
relaxation processes we only have taken into account p-p 
Coulomb interaction. This assumption should be justified 
and for that purpose we are looking for some alternative 
processes which could support the relaxation of the proton 
distribution function. 

Since the proton distribution function resulting in the 
interface is a consequence of charge-exchange induced in- 
teraction between the distributions of protons and H- 
atoms, the proton distribution function f p approximately 
represents a superposition of two shifted Maxwellians with 
different partial densities uq and n\, the latter being due 
to ionized H-atoms. Under such circumstances it is sug- 
gestive to first investigate whether excitation of electron 
plasma waves by the positive slope of the composite dis- 
tribution function may be an important contribution to a 
relaxation of such non-equilibrium distributions towards 
a best-adapted mono-Maxwellian. 

Taking the approach presented by Davidson et al. 
(1970) and going to the limit of full equilibration of elec- 
tron thermal and differential proton bulk energies, one ob- 
tains a friction force i? e ,i per unit volume acting between 
protons and electrons given by (see Fahr & Neutsch 1983): 



Re,i = ^TmaxWfeoy J^T (63) 

where J = 39.9 is found as a constant, m e and T e are the 
electron mass and temperature, respectively, 7 max is the 
maximum growth rate of the excited plasma waves (i.e., 
most unstable electrostatic mode) given by: 

a/3 / m e 

^^ife (64) 

where tu p is the electron plasma frequency, and Wko is 
the total wave power which according to Davidson et al. 
(1970) in the saturation phase attains the following value: 

W k0 ~ n e m e (u H - u p ) 2 (65) 

where uh and u p are the z-componcnts of the local bulk 
velocities of the H-atoms and the protons, respectively. 
The above mentioned friction force R e j operates with the 
tendency to reduce the differential bulk kinetic energy in 
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the proton distribution function. Thus one can derive a 
typical relaxation time r e ^ with the following relation: 



ni(\m p {u H - u p ) 2 ) 



(66) 



Inserting values like no = 0.1 cm -3 , m/no =0.1, uh = 
25 km/s, Up = 10 km/s, and KT e — jm e (uH — u p ) 2 one 
gets the result: r e ^ ~ 2000 s. 

From the above estimate one should expect to have 
the relaxation of the non-equilibrium proton distribution 
occurring via electron plasma wave excitations within as- 
tonishingly short times r e> i of the order of one or a few 
hours. 

One should, however, check the validity of the assump- 
tions made in the above derivation of the time r e ,i. As it 
turns out, the results of Davidson et al. (1970) can in 
principle be adopted only if the two Maxwellian peaks in 
the superposed distribution function are clearly separated 
from each other and do not essentially overlap by their 
Maxwellian wings. When looking at the conditions which 
in fact are prevailing at least on the interstellar side of the 
heliospheric interface it becomes fairly evident that this 
condition is by far not fulfilled, since the proton thermal 
velocity is of the order of 20 km/s and the separation of 
the bulks only amounts to perhaps 15 km/s. 

This conclusion can also be supported by looking 
into the Landau damping rates or the plasma oscillation 
growth rates 7 P (w, k) which are connected to the deriva- 
tives with respect to v z of the normalized distribution 
function f£ z — J J fpdv x dv y by the following relation: 



7 P (w,fc) 



p 

irk 2 



dv z 



v z —uj I k 



d fp. 



dv. 



(67) 



where the normalized distribution function /* may ap- 
proximately be represented by: 



r X / \ r 21 1 Tl\ I Tp 

fp,zK w z) ^ ^exp [-w z \ + —j=— -f- 



1/2 



(68) 



x exp 



-^(W H -W p -w z ) 

J-H 



wt 
W 2 



Here the following denotations were used: 
m p v 2 /(2KTp); W 2 = m p u 2 H /(2KT p ); 
m p u 2 /(2KTp); T p and Th are the temperatures of the 
bulk of the protons and of the H-atoms, respectively. Now 
from the above relations one finds: 



dw z 



-2w, 



■ exp|— w z 



2 m 



T 



3/2 



(69) 



x (Wh — W p — w z ) exp 



Tp 
T) 



H 



w, 



The most positive contribution comes from the second 
hump in the distribution function to fp Z (w z ) and leads to 
positive growth rates evidently coming from regions of w z 
where the following relation is valid: Wh — W p — w z ~ 1, 



which then leads to: 
i 1 

2 m 



dw z 



W v 



■exp[-(W H -W p -l) 2 



H 



3/2 



exp 



T 

T] 



H 



(70) 



Taking now typical values for the above quantities as 
they prevail in the LISM interface we can decide whether 
positive growth rates are possible. In view of the val- 
ues we display in Figures and [21 of Section 1 we se- 
lect the following typical values: Wh — 3; W p ~ 0.3; 
T p /T H ^ 12000/8000 = 1, and thus: 



dw z 



exp[-(1.7) 2 ] + A^ ( 1.5)le-i (71) 



exp[-(1.7) z ] ( — (1.5) 3/2 exp[+0.2]-l 



n 



meaning that only under unlikely conditions, not met in 
our calculations above and characterized by: 



ni 
no 



> 



1 



1.5 3 / 2 



exp(-0.2) = 0.44 



(72) 



a relaxation of the distribution function f* z by excitation 
of electron plasma waves is likely to occur. 

It thus remains to inspect relaxation processes con- 
nected to the nonlinear interaction of protons with MHD 
turbulences, amongst which pitch angle diffusion processes 
seem to be the most effective ones. These, however, de- 
pend on the local magnetic fields and the existing MHD 
turbulence levels. Of course, nothing is known at present 
about LISM magnetic fields and turbulences. This is why 
no good estimate of pitch-angle diffusion processes is pos- 
sible here as is otherwise feasible on the basis of an expres- 
sion already derived by Hasselmann & Wibberenz (1970) 
in which the typical time period for pitch angle diffusion 
is given by: 



V 



+1 



{l-u 2 fdn 



D„ 



(73) 



where Am is the mean free path with respect to pitch- 
angle scattering, \i is the pitch angle cosine, and -D w is 
the Fokker-Planck pitch angle diffusion coefficient. The 
latter, however, is a complicated function of particle ve- 
locity, local magnetic field magnitudes, Alfven velocities, 
and turbulence levels. This coefficient D^^ then leads to 
velocity-dependent mean free paths An(i;) and to corre- 
sponding scattering periods T|| (v) (as given by Chalov & 
Fahr 1999a, b for keV-energetic ions). The needed input in- 
formation to evaluate the above expression unfortunately 
are completely lacking for the LISM interface region and 
furthermore the newly injected ions from the H-atom pop- 
ulation in this region only have energies of about 10 eV. 
Thus we cannot estimate tii in this way. 

Instead we estimate the relaxation times connected 
with pitch angle scattering on the basis of the reasonable 
expectation that the locally appearing MHD turbulence 
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levels are not primarily due to convected pre-existing tur- 
bulences, which have not been quantitatively described 
in the literature, but due to wave-driving processes of sec- 
ondary protons injected via charge exchange processes into 
an unstable mode of the proton distribution function. Let 
us assume that the newly created protons produced by 
charge exchange are populated in a region of velocity space 
where they represent some free kinetic energy by which 
they are able to drive some MHD wave power (see, e.g., 
Huddleston & Johnstone 1992, Williams & Zank 1994, 
Fahr & Chashci 2002). Let us then assume that this free 
energy is pumped into turbulent wave power essentially at 
an injection wave number fcj. The following considerations 
are of course only relevant if an interstellar magnetic field 
is present. For an unmagnetized interstellar medium what 
follows here would thus be irrelevant. 

Adopting that near the symmetry axis (stagnation 
line!) a background magnetic field Blism exists which 
is oriented parallel to the stagnation line (i.e., according 
to an MHD interaction model with Blism \\ V lism as 
treated, e.g., by Baranov & Zaitsev 1995, Pogorelov & 
Matsuda 1998, or Ratkiewicz & McKenzie 2003), one then 
can assume that Alfven waves are excited by this free en- 
ergy which then propagate upstream and downstream of 
the plasma flow. The free energy of newly injected pro- 
tons stemming from the bulk of the H-atom distribution 
function can then be assumed to be pumped into the wave 
field at the injection wave number fulfilling the resonance 
conditions and given in the bulk plasma system by: 



=F h ~ n/ (v A ± (UH - Up)) 



Since hydromagnetic Alfven waves have frequencies 
u> = va k < fl the cyclotron resonance condition can only 
be met by ions with positive/negative parallel velocities 
resonating with waves with negative/positive wave vec- 
tors. In the outer interface plasma newly injected ions with 
uh > Up have positive parallel velocities in the plasma 
frame and thus can only resonate with waves propagating 
in the direction opposite to them with wave vector magni- 
tudes h ~ Q / {v a+{uh—u p )) . When resonating with these 
waves, energy can be exchanged between such ions and the 
waves up to an equilibrium situation. This injected wave 
power then cascades from fcj to larger wave numbers de- 
scribed by a so-called wave-wave diffusion coefficient Dk k 
(see Zhou & Matthaeus 1990) until the power arrives at a 
wave number kdi S ~£1/va where it is resonantly absorbed 
by the bulk of the protons. 

We now determine the time Tkk needed to cascade the 
injected power to the dissipation wave scale - in order 
thereby to derive the time needed to transfer energy from 
the secondary proton hump in the proton distribution 
function to the primary proton bulk. Following Chashei 
et al. (2003) we obtain: 



i~kk 



(Afc) 2 _ (hi - k dls y 



VA + (uh—U p ) 



VA ) 



AD 



kk 



AD 



u, 



4D 



kk 



(75) 



According to Zhou & Matthaeus (1990), we can represent 
the k-space diffusion coefficient by the following expres- 



7/2 4nWj_ 



D kk =Cv A k'r\l-^- (76) 

£>LISM 

where Wi denotes the spectral power at wave number ki, 
and C ~ 0.1 is a constant. Assuming that wave energy is 
injected at the same rate as it cascades down to kdis, we 
can determine Wi by the following relation: 

Afc 2 

WiAfc ~ PexAe^kk = /3exAe;— — (77) 

4L>kk 

where (3 ex = n v nn<j ex V re i is the mean charge ex- 
change rate between H-atoms and protons, and Ae,; = 
ef [\rn p (uh — Up) 2 ^ is the free energy pumped from the 

newly injected protons into the wave field. With Eqs. I|75|l 
through (|77|l one thus finds: 



W, 



3/2 



Blism a A » , v fpp~f3exAe i Ak 

/3 ex AsiAk = - ( 78 ) 



2Cv A kJ /2 V^ 



Ck 



With the above result we finally find the diffusion time 
T kk as given by: 



i~kk 



W,Ak Afc / ^fpp-pexAeiAk 



2/3 







4Cfc 

5/3 



(79) 



(74) 2C 2 / 3 (n H VexV re l c f (u H - u p ) 2 kfj 



1/3 



Let us now assume for estimation purposes that V re i 
can be reasonably well approximated by V re \ ~ (uh — Up). 
Then we arrive at the following expression: 



Tkk 



i 

[uh ~ Up) 



(80) 



2/3 



2CV^vT (riff a ex e f kff 3 

_ (UH ~ Upf /3 {VA + U H - Up) 

~ 2e 1 / 3 C^v 5 fn^(n H a ex ) 1/3 

For further evaluation of this expression we transform it 
into the following form: 



Tkk 



(X H -X p ) 2/3 (l+X H -X p ) 2/3 
2e 1 / 3 C 2 W3 (v A n H a ex ) 1/3 



(81) 



where the quantities Xh and X p denote the corresponding 
bulk velocities normalized by va- 

It remains to find some adequate value for the frac- 
tion Cf of energy that is transferred as free energy to the 
wave fields. This value we intend to derive in analogy to 
a similar derivation given by Chashei et al. (2003) for the 
free energy of newly injected pick-up ions in the inner he- 
liosphere. The following estimation cannot be carried out 
prior to at least some qualitative view of the configuration 
and magnitude of the magnetic field Blism in the outer 
interface. 
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Hence we assume here that close to the stagnation line 
on the interstellar side of the interface we have a field 
-Blism which is parallel to the plasma flow, i.e., parallel 
to u p (see, e.g., the model presented by Baranov & Zaitsev 
1995). Then newly ionized H-atoms are implanted in the 
proton distribution function roughly at a velocity Am = 
\uh — u p\ relative to the proton bulk flow. Injection of a 
newly created proton from the H-atom bulk on the average 
will eventually lead to the following relative velocity g with 
respect to the proton bulk: 



g~yjAv? + (? H 



(82) 



where ch is the thermal velocity of the H-atoms. 

Now we consider the two frames of Alfven waves mov- 
ing with the velocity ±v A relative to the proton bulk frame 
in up-field and down-field directions. With respect to these 
wave frames the injected protons have relative velocities 
given by: 



9A± = y(v A ± Auf + c 2 H = x /v A ± 2v A Au + g 2 (83) 

and by resonance with the up-field waves will rapidly 
be pitch-angle scattered onto the accessible spherical 
shell that is associated with the upgoing wave-frame (see 
Huddleston & Johnstone 1992). With respect to the pro- 
ton bulk frame they thereby lose an energy Ae/, taken 
as the pitch-angle average, given by (see Fahr & Chashci 
2002, Chashei et al. 2003): 



(84) 



cos 



0+ 



9A+ 



(86) 



After evaluation of the integral in the above expressions 
one then finds: 



(g% {■&)) = 2v\-2v A Au + g 2 



v A c 2 H 



y/v A + g 2 - 2v A Au 



(87) 



which leads to: 

Ae f+ ~ \m p [g 2 -(gl(d))} = 



(88) 



--m p g 2 ( 2z\ - 2z A z v 



z A z H 



y/l + z A - 2z A z v 



where the following notations have been used: z A — v A /g; 
z H — Cn/g; and zjj = Au/g. If we now put into the 
expressions above some concrete numbers, wc find: 



v A = v A0 (B LIS m/B ) y/n /n L isM 



(89) 



where the quantities with index "0" are reference quanti- 
ties and, taking these reference values for the solar wind at 
the orbit of the earth, i.e., v A o = 40 km/s, B = 5 Gamma 
and no = 5 cm" 3 , we finally obtain: 



v A = wao(3/50)V5/0.1 ~ 17 km/s 



(90) 



Then taking U H ~ 20 km/s, U p ~ 10 km/s, and c H ~ 
15 km/s, we obtain: 

g ~ ^AU 2 + cjj = V10 2 + 15 2 km/s = V325 km/s 
= 18 km/s; z A = v A /g = 17/18 = 0.94, z v = AU/g = 
0.55; z H = c H /g = 0.83; X H = U H /v A = 20/17 = 1.18 ; 
X p = U p /v A = 10/17 = 0.58 and thus wc find: 



Ae /± 



1 



m p g 2 (2 • 0.94 2 - 2 • 0.94 • 0.55- 



(91) 



0.94-0.83 2 



VI + 0.94 2 - 2 -0.94- 0.55 



= 0.55 



which due to the fact that newly created protons only 
resonate with waves propagating in the up-field direction 
finally yields e/ + ~ 0.55. Hence one obtains: 



Tkk 



(1.18 - 0.58) 2/3 (1 + 1.18 - 0.58) 



2/3 



2 • e^C 2 /^ 2 /* (v A n H a ex ) 



1/3 



leading to: 



Tkk 



(0.6) 2/3 (1.6) 2/3 



Ae f ~ X -m v [g 2 - (g\ (d))] = e f Qm p g 2 ^j 

where the angle-averaged velocity (<7+(#)) of the pitch- 
angle scattered new protons is given by: 
l 

/ (v A + g A+ - 2v A g A ± cos i?) d cos t? 
(g\ W ) = p (85) 

1 - COS #Sax 

Here the maximum permitted inclination angles are de- 
fined by: 

v A + Au 



2-0.55 1 / 3 C 2 / 3 fl 2 / 3 (n H v A a ex ) 



1/3 



Putting in numbers like: C — 0.1; nn = 0.1; <7 e 



(92) 



(93) 



2.5 



9.5 • I0- 2 (B LISM /Bo) = 5.7 • 10~ 3 s we 



then obtain: 



Tkk 



0.71 • 1.36 



2 • 0.82 • 0.22 • 0.033 (1.6 • 10 6 • 2.5 • 10- 16 ) 1/3 
10 4 s 



(94) 



This now means that the isotropization and redistribu- 
tion of newly incorporated protons from the primary ring 
distribution to the associated spherical shell distribution 
will occur within typical time periods of a few hours. 
This, however, does not mean that a full relaxation to- 
wards an associated equilibrium distribution is likely to 
occur within a similar time period. The resulting distri- 
bution function / is still far from fulfilling the require- 
ment 6 (J / In / d 3 v) /St = 0. One other point, however, 
becomes very clear when comparing the period Tkk with 
the period t p yielding T p /Tkk ~ 10 4 , namely: the process of 
wave driving operates much faster than the proton-proton 
relaxation, so that at least the fraction £/ of the kinetic 
energy of newly implanted ions is pumped into wave tur- 
bulence. 

5. Conclusions 

In purely hydrodynamic multi-fluid interaction codes de- 
scribing the interaction of the partially ionized interstel- 
lar medium with the solar wind plasma, both the proton 
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fluid and the H-atom fluid are generally represented by the 
moments of distribution functions which are taken to be 
bulk- flow shifted Maxwellians. We have shown for the re- 
gion close to the stagnation line, outside of the hcliopause 
but inside of the outer bowshock, that this assumption is 
substantially violated for the H-atom flow, while it is only 
mildly violated for the proton flow. This result, however, 
is due to the values of interstellar parameters adopted in 
our study, where we have taken the LISM H-atom density 
= LISM proton density = 0.1 cm~ 3 . In that case the p-p 
Coulomb relaxation processes are quite effective in keep- 
ing the resulting proton distribution function close to a 
shifted Maxwellian. 

For different LISM parameter combinations, e.g., for 
lower LISM proton densities, or at different places of the 
heliospheric interface, the assumption of having the proton 
distribution function close to a shifted Maxwellian may be 
violated. This is because the Coulomb relaxation period 
t p ~ Tp^ 2 /n p which means that for lower proton densi- 
ties and higher proton temperatures this relaxation period 
may easily increase to values t p >> r ex and thus induce 
substantial deviations of the resulting proton distribution 
function from a shifted Maxwellian. For instance, inside 
the hcliopause, where very low proton densities and high 
proton temperatures prevail, substantial deviations of the 
proton distribution from shifted Maxwellians should be 
expected unless alternative relaxation processes operate, 
as discussed in Section 4. 
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Fig. 4. Parameters relevant to compute the proton relaxation time, discussed in Ea. 1531 through IfTTl as functions of 
the heliocentric distance between the bowsliock and the heliopause, based on the parameters shown in Fig. |3 
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Fig. 5. Hydrogen and proton distribution functions: no proton relaxation case. Light gray lines correspond to protons 
and dark gray to H-atoms. Shown at selected distances from the Sun (in the region between the bow shock and the 
heliopause) are the unperturbed Maxwellian distribution functions, drawn with broken lines, and the net distribution 
functions f p -n — fa,p;H + /i,p;£f ; drawn with the solid lines and expressed in atoms cm -6 s -3 . 
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Fig. 6. Ratios of the perturbation functions fi, P -H to the corresponding unperturbed distribution functions fo,p,H' the 
/i//o ratios, in the case of no proton relaxation operating, drawn for the H-atoms in dark gray and for the protons in 
light gray. 
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Fig. 7. Hydrogen and proton distribution functions: proton Coulomb relaxation case. Light gray lines correspond to 
protons and dark gray to H-atoms. Shown at selected distances from the Sun (in the region between the bow shock 
and the heliopause) are the unperturbed Maxwellian distribution functions, drawn with broken lines, and the net 
distribution functions f p -H = fo,p;H + fi,p;H, drawn with the solid lines and expressed in atoms cm -6 s~ 3 . 
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Fig. 8. Ratios of the proton perturbation function to the proton unperturbed distribution function: /i//o for the 
proton Coulomb relaxation operating. 



